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Abstract: We study Q-ball type solitons in arbitrary spatial dimensions in the 
setting recently described by Kusenko, where the scalar field potential has a flat di- 
rection which rises much slower than 2 . We find that the general formula for energy 
as a function of the charge is, Ed ~ Q d d ^ d+1 ^ in spatial dimension d. We show that the 
Hamiltonian governing the stability analysis of certain Q-wall configurations, which 

i__> 

are one dimensional Q-ball solutions extended to planar, wall-like configurations in 
three dimensions, can be related to supersymmetric quantum mechanics. Q-wall 
and Q-string (the corresponding extensions of 2 dimensional Q-balls in 3 spatial di- 
mensions) configurations are seen to be unstable, and will tend to bead and form 
planar or linear arrays of 3 dimensional Q-balls. The lifetime of these wall-like and 
string-like configurations is, however, arbitrarily large and hence they could be of 
relevance to cosmological density fluctuations and structure formation in the early 
Universe. 
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1. Introduction 

Non-topological solitons were first described in the mid-70's Jl[] in the context of 
scalar field models with two fields, symmetry breaking and a conserved charge. They 
reappeared in a somewhat different guise in the mid 80's in theories without spon- 
taneous symmetry breaking as considered by Coleman and the name Q-balls was 
invented. In the mid 90's they have re-surfaced in connection with models of scalar 
fields where the potential has flat or quasi-flat directions, which is generic in models 
of supersymmetry breaking, as described by Kusenko @. 

The original non-topological solitons and Q-balls had an energy which behaved 

as 

E~fiQ (1.1) 
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for large conserved charge Q. If /i is less than the mass of a perturbative, charged 
excitation of unit charge, then the Q-balls are clearly stable against decay into per- 
turbative charged particles. Indeed, Coleman defined Q balls as the minimal 
energy classical configurations of fixed charge Q and was able to show their absolute, 
classical stability. However, it was soon realized that linear terms induced via quan- 
tum corrections in the effective potential would lead to the decay of these Q-balls 

I- 

The recently-described Q-balls do not suffer this fate. They have an intriguing 
energy-charge dispersion, E ~ Qw 4 ) for large Q in 3 spatial dimensions. Hence 
these Q-balls are arbitrarily tightly bound compared to the energy of an aggregate 
of perturbative quanta of equal charge. We note that such an asymptotic energy- 
charge dispersion with the non-integer power 3/4 also arises in another context. For 
solitons in models with non-trivial topological Hopf windings 0, there is an energy 
bound, which is not attained, but behaves as Q^ 3 ^\ It would be interesting to find 
a model where a "Hopf soliton", perhaps of the Q-ball, type can exactly saturate 
such a bound. Indeed, we will find below exactly solvable Q-balls which do saturate 
a similar bound in 1 spatial dimension. 

In this paper we investigate extended objects of the Q-ball type. These are the 
solitons of one- and two-dimensional models embedded in three dimensions, forming 
Q-walls and Q-strings, respectively. Similar structures have been studied recently 
in |IT[ both analytically and using numerical simulations. The cases discussed 
concern Q-balls which have an energy-charge dispersion that grows linearly with the 
charge, E ~ [iQ. It is found that collisions of such Q-balls can give rise to Q-string 
like states. In JT2] stable Q-ring solitons are described. Here stability is imbued via 



an additional topological twisting of the phase as one goes around the ring. We find 
our study is complementary to these studies. 

In the next section, we study Q-solitons in arbitrary spatial dimension d. The 
energetics of these objects are interesting: the energy Ed of a soliton in d dimensions 
depends on its charge Qd as Ed ~ Q d J . The lower the dimension, the more 
energetically favourable the Q-soliton over a conventional configuration of charge 
Qd- At first glance, this would seem to indicate that extended objects in three 
dimension would be energetically preferred over Q-balls. However, as we shall see, 
this turns out not to be the case. 

In Section 3, we consider embedded lower-dimensional Q-solitons in three di- 
mensions. The energy E 2 and charge Q2 of a two-dimensional soliton, for instance, 
are the energy and charge per unit length of the corresponding Q-string. The energy 
of a Q-string of length L and of fixed three-dimensional charge varies as L 1 / 3 . Thus 
the tension varies as L _2//3 , reducing rapidly with length. However, the energy per 
unit length also varies as L -2 / 3 , so the natural time scale for oscillations of Q-strings 
is equivalent to that of topological strings (whose energy per unit length and tension 
are independent of L). 
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In Section 4 we study the stability of Q-strings or walls against "beading", i.e., 
the formation along the length of the string (or surface of the wall) of clumps of higher 
charge (and energy) per unit length (or surface) which would ultimately result in the 
formation of a one-dimensional chain (or two-dimensional sheet) of non-interacting 
balls. 

We study the stability of walls, in particular, since there the analysis is much 
simpler, essentially because the one-dimensional soliton can be found analytically 
for certain potentials. In studying small oscillations about such a solution, the 
eigenvalue equation determining the frequency of these oscillations can be related 
to a Schrodinger equation in supersymmetric quantum mechanics. Applying known 
results from that theatre enables us to identify a single unstable mode for the fluc- 
tuation, two zero-modes (one of which corresponds to translation of the wall) and a 
spectrum of stable modes. The unstable mode can have an arbitrarily large charac- 
teristic time, so the wall can be extremely long-lived. We expect the same behaviour 
for strings, but have not yet done the analysis in that case. 



2. Q-balls in any dimension 

Consider a charged scalar field in d spatial dimensions, the Lagrangian governing its 
behaviour is 

C = j d d X (^d^ a d^ a - V(\<t>\^J (2.1) 

where a = 1,2, |0| = a/0? + 02 an d we use a simple SO(2) global gauge invariance 
whose conserved charge corresponds to Q. The conserved charge is normalized as 



Q 



J d d xe ab (t)adth- (2.2) 



The scalar potential, inspired by supersymmetric models, satisfies 

where A is a constant which is determined by the scale of supersymmetry breaking. 
The minimum of the potential occurs at |0| =0 and the charge symmetry is unbroken. 
Q-ball solutions are found with the ansatz 

1 + i0 2 = e-V(O (2.4) 

where ui is a constant, r is the radial coordinate and <p(r) is real. We immediately 
obtain 

Q = u J d d xcp 2 (r). (2.5) 
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The profile taken for <p(r) satisfies the asymptotic conditions 

■*) = {+> rSS0 (2.6) 

where 0o is taken large enough so that the asymptotic behaviour of the potential 
is valid. Hence the Q-ball interpolates from large non-vacuum values around 4>q in 
its center to ip = 0, the perturbative vacuum, at infinity. The existence of stable 
solutions to the equations of motion for <p(r) in 3 spatial dimensions is the content 
of the previous papers || |] and is easiest seen by the classic continuity arguments 
of Coleman 0. Here we re-do the analysis in arbitrary spatial dimensions d. 
The equations of motion for (p(r) are 

VV(r) - V'{ip{r)) + u 2 (p(r) = (2.7) 

which, given the form of the potential, become 

V 2 ip(r) + uj 2 ip(r) = r^O 

VVO) - m 2 (p(r) = r^oo (2.8) 

using the assumed asymptotics of the potential. The solutions of these equations are 
immediate, but dimension dependent. For small r we get spherical Bessel functions 
of order in 3 dimensions, ordinary Bessel functions of order in 2 dimensions, and 
trigonometric functions in 1 dimension, with argument uor. For large r the argument 
simply changes to imr, which is imaginary. 

2.1 d = 3 

The asymptotic solutions are 

„ W = (2.9) 
[ ip k {mr) r oo 

where 0o is a constant that is fixed by the total charge of the solution, ip is a 
dependent constant and joiur) and k (mr) are the spherical Bessel function of order 
of real and imaginary argument respectively. The asymptotic solutions interpolate 
from one to the other in a relatively small region, and it is difficult to obtain analytic 
solutions, except in the 1 dimensional case. However, the continuity argument of 
Coleman easily establishes the existence of a solution. 

The idea is to first observe that the equation of motion |2.7j can be interpreted to 
describe the classical motion of a "particle" in one dimension of position ip, in the 
potential — V(ip) + \io 2 ^p 2 and in the presence of a specific frictional term that is time 
dependent coming from the Laplacian, and where r plays the role of time. Indeed 
the equation of motion is equivalently written as 

(p(r) - VUpir)) + wV(r) = --<p(r) (2.10) 

r 
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where (p{r) = d r (p(r). The potential is negative for small (p, and is approximately a 
downward opening parabola of curvature — m 2 , whose maximum passes through the 
origin at <p(r) = 0. It turns over, after achieving a minimum (negative) value, into 
an upward opening parabola of curvature u 2 for large enough values of (p(r). Clearly, 
eventually the potential becomes positive. Now if we start with an initial value 0o for 
which V(4>o) is negative, the "particle" will never rise up high enough to reach the 
origin, even in infinite time, and will perform damped oscillations till it settles down 
to the global minimum value of —V(<p) + \<J 2 ^> 2 . On the other hand if we start with 
a large initial value of <po, then since the "motion" is governed arbitrarily well by the 
spherical Bessel function jo(ur), the "particle" will attain a given value close to the 
point where the full potential — V(ip) + \uj 2 (p 2 vanishes, in a fixed time, independent 
of the initial value 0o which is of the order of 1/u. The velocity at this point will be 
proportional to 0o- Hence the subsequent evolution can be arranged to commence 
with an arbitrarily high initial velocity by choosing 0o larger and larger. It is clear 
then, that there exists an initial value for (f) so that the "particle" will arrive at this 
point with sufficient velocity that the subsequent motion will rise up and go over 
the maximum at (p — in finite time. Hence, by continuity, there exists an initial 
value for 0o which will exactly achieve the maximum as time goes to infinity, i.e., 
for r — > oo. Such a configuration is the solution and it will give the appropriate 
interpolation between the two asymptotic solutions pT9| . Clearly the value of ip will 
be governed by the value of 0o which in turn will be fixed by the value uj. Hence the 
total charge Q can be thought of as a function of u alone, or conversely all of the 
parameters depend only on the total charge Q. 

The energy of such a configuration can be approximated as 



where we neglect the contributions from the surface and the exterior of the Q-ball, 
and where a and (3 are positive constants. This expression is of course reasonable 
only for Q-balls with large charge Q. Using 



E 




(2.11) 




(2.12) 



where 7 is another positive constant, we get 



E 



oiQuj Aj3 



(2.13) 




Imposing dE/duj = yields 



aQ 



3A(3 



(2.14) 
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which gives 



u ~ Q- 1/A 

00 ~ g 1/4 

£ ~ Q 3/4 . (2.15) 



2.2 d = 2 



Here the solution for small r is 0o Jo(^ r ) where Jq{ujt) is the ordinary Bessel function 
of order 0, and for large r the solution is ipoK (mr), where K (mr) is the Bessel 
function of order zero of imaginary argument. The energy is then approximately 
given by 



E 2 « / d 2 x (u; 2 <p 2 (r) + \Vif(r)\ 2 + A) 

Jr<l/w 



= ^ + ^ (2.16) 
which yields upon imposing dE 2 /du = 

00 ~ <?2 /3 

^2 ~ QT- (2.17) 

Here i? 2 and Q 2 are the appropriate 2 dimensional energy and charge respectively, 
which in terms of 3 dimensional energy and charge are the energy per unit length 
and the charge per unit length. 

2.3(1 = 1 

Here the equation of motion can be reduced to quadratures, and it can be integrated 
exactly in many cases. It is more convenient to use the normal spatial coordinate 
x G [— oo, oo] rather than the radial coordinate r. The equation of motion for <p(x) 
where 4>\(x) + ifaix) = e lult (p(x) is 

which yields 



/ 

Jd> 



= / dx. (2.19) 

lc/,0 ±y/V((f) ~ Ld 2 (f 2 J xo 

Here x gives the position of the Q-ball and the ± is chosen depending on whether 
x > Xo or x < xq. The square root should not become imaginary thus we must have 
that V((p) > u 2 (p 2 . Since for large values of (p the potential becomes flat and the 
required inequality is violated, there is a fixed maximum value for O given by the 
condition V((j>o) = u; 2 0q. There is no "friction" term coming from the Lapalcian in 1 
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dimension, hence conservation of energy dictates that this value of 0o is the starting 
point for (po(x). Small to implies large values of 0o which gives approximately that 



;,2 



m A/u 2 that is 0o ~ 1/uJ. The quadrature [2.19| is not particularly useful in the case 



of a general potential V(ip(x)). The solution for x ~ x is ip(x) = 0o cos(a;(a; — x )) 
while for |x — | >> n /2uj is y?(=bc) = ifoe Tmx . The charge is given by 



xq+(-k/2ui) 



i7T 



Qi = io / <iry2 (x) ~ uj(f) Q / dx cos (oo(x — xo)) = — — (2.20) 

J Jx -(tt/2lu) 2 

while the energy is 

E 1 = [ dx{u 2 <p 2 (x) + (DM*)? + V(<p(x))) 
dx(2u 2 ip 2 + V{ip{x)) - ip{x)V'{ip{x))) 



Imposing 

yields 

implying 

This gives, 



dx(2u 2 ip 2 + A) 

\x— xo\<(tt/2lu) 

2cu 2 0q( / dx cos 2 (lu(x — x )) ) H 

2^| + Att_ (221) 
2 

do; 



(2.22) 



i j_A = ao 1 _A =0 



= ^. (2.24) 



2Qi 



£ 2 ~ QT- (2.25) 
2.4 Arbitrary spatial dimension 

Here we record the corresponding formulae in arbitrary spatial dimension d. The 
Qrf-ball solution will be to the appropriate generalization of the free spherical wave 
in its centre with argument ur while from its edge it will evolve into the exponentially 
decaying type solution of mass m outside, i.e., the same function but of imaginary 
argument imr. One easily finds 



Qc 



iO d Co^- 1 ) 
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E d « a^J 2 ^ + M 

a d Q d uj( d -V , 2 _ d) (3 d A 

— Ld x H -T- 

= *Q*L + M (2 . 26) 
Id oj d 

where a d , f3 d and 7^ are the corresponding constants in d dimensions. Imposing 
dE d /du = yields in a straightforward manner 

E d ~ Qj /(d+1) . (2.27) 
3. Energetics of Q-strings and Q-walls 

Suppose initial conditions have created a large Q-wall or Q-string in 3 dimensions. 
We will observe here that such structures have interesting energetics. We address 
the question of stability of such objects in the next section. 

3.1 Q-strings 

We define Q-strings as linear configurations in 3 dimensions whose cross section 
corresponds to a 2 dimensional Q 2 -ball, with a fixed charge per unit length Q 2 . Such 
a configuration will simply terminate at two end points or close on itself yielding a 
closed Q-string configuration. We will consider for convenience the second situation, 
of a closed Q-string of length L. We will assume that L is sufficiently long so that 
gradient energies due to the gradual bending that is necessary to close the Q-string 
are negligible. The relevant conserved charge of a closed Q-string then is given by 

Q = Q 2 L (3.1) 

while the energy is given by 

E = E 2 L. (3.2) 
Replacing for E 2 ~ Q% 3 = (Q/L) 2 ^ 3 we get 

E ~ Q 2 2 /3 L = Q 2/3 L 1/3 . (3.3) 

Hence a Q-string of total charge Q satisfies an energy-charge dispersion relation 

^-string ~ Q 2/3 - (3.4) 

This does not mean that the Q-string is even more tightly bound than a Q-ball with 
the same charge, which satisfies E ~ Q 3 ^. Overall scaling of the Q-string always 
lowers its energy. The energy-length dispersion is 

Eq —string ~ 

l}l 3 (3.5) 
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so that as the length decreases, so does the energy. Indeed the Q-string will shrink 
and reduce its energy until its length is comparable to its width, at which point the 
description in terms of a Q-string is no longer valid. 

The energy-length dispersion is however, remarkable. The restoring force which 
tends to cause the Q-string to shrink 

dE IQ 2 / 3 

String - ~ - 3^3 (3.6) 

disappears as L gets very large. This is in contradistinction to the case of a topologi- 
cal string, such as a cosmic string or vortex loop, where the energy-length dispersion 
is -Etop.-string = fiL, where \i is the mass per unit length, giving a constant restoring 
force. Nevertheless, the time scale for shrinking of a Q-string is not markedly differ- 
ent, for a large Q-string the mass per unit length also decreases accordingly with its 
length, compensating for the reduced restoring force. 

The time scale for shrinking of a topological string is governed by the effective 
Lagrangian 

1 L 2 LL 2 

Aop.-string = 2 ^4^2 ~ ^ = ^~g^2 ~ ^) (3-7) 

while for a Q-string we have 

L -V g ML^ ( aQ23 )L-a0 2/3 ( Ll/3p L^) (3 8) 

^Q-string - 2^^W } 4^ ~ (~tfJT } ~ aQ2 ( S7T 2 ~ ^ ' 

We see that neither /x nor aQ^ 3 play a role in the dynamics. Conservation of energy 
for the topological string gives 



while for the Q-string 



8tt 2 
L^L 2 



LL 2 

+L = L (3.9) 



g7r2 -L^ = L]J\ (3.10) 

where Lq is the initial length in each case. With the change of variables L = v 2 for the 
topological string and L = u 6 for the Q-string the equations become easily integrable, 
yielding the time for contraction T top _ string = L /A^ and TQ_ str i ng = 45L /32\/2. 

We will see in the next section that the Q-string and the Q-wall are both unstable 
to beading into arrays of smaller Q-balls. If the Q-string were to spontaneously decay 
into a linear array of N, Q-balls each of charge q, we have 

Q = Nq (3.11) 

from conservation of charge, while the energy of such a final state is 

E = (3q 3 ^N = (3Q 3/4 N 1/4 . (3.12) 
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Evidently such a state has higher energy than the absolute minimum of energy in 
the sector of charge Q, given by the Q-ball of charge Q, (N — 1), which has energy 
E = f3Q 3 ^. But it will be an intermediate state with less energy than the Q-string 
if (adopting the notation of Section 2.4 to avoid confusion) 

(3 3 Q 3/4 N 1/A < p 2 Q 2 / 3 L 1 / 3 . (3.13) 

This implies 

N < <£)<|£ (3,4, 
which gives N/L < (/V/3 3 ) 4 Q 2 " 1/3 . Asthe 

size of each smaller Q-ball is approximately 
R ~ L/N thus R > (fa / (32) A Q^ 3 , i.e., the Q-string will be stable against decay into 
a chain of N, Q-balls if the size of each Q-ball is smaller than Q l J 3 ~ 1/cj, which is 
just the width of the original Q-string. 

3.2 Q- walls 

We can repeat the above analysis for the case of a closed Q-wall of total 3 dimensional 
charge Q and with a surface area L 2 . Then 

Q = L 2 Q 1 (3.15) 

and 

£ Q -waii = L 2 E, = L 2 f3 1 Q\ /2 = (3iL 2 (Q/ L 2 f 1/2) = ^Q^L. (3.16) 
This is to be compared with the energy of a topological domain wall, 

-^top. -wall = ^L 2 (3-17) 

where \x is the mass per unit area. For a spherical surface the effective Lagrangian 
describing the dynamics of a spherical topological domain wall is 

1 L 2 

Aop.-waii = -zfJ-L fiL (3.18) 

Z 47T 

while for a spherical Q-wall we get 

/Vwaii = l^Q 1/2L ^ ~ PiQ 1/2 L. (3.19) 
These yield the conservation law for the topological domain wall 

-^—L 2 L 2 + L 2 = L 2 (3.20) 

07T 



while for the Q-wall 



^-LL 2 + L = L , (3.21) 

07T 
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where L is the initial value of L. These are easily integrable and yield the time for 
contraction to be 

^top.-wall = r— (3.22) 

2v 27r 

and 

Tg-waii = \J~~^~^ (3.23) 

for the topological domain wall and for the Q-wall respectively. Again we see that 
the time for contraction is about the same in the two cases. 

If the Q-wall spontaneously beads up into a locally planar array of N smaller 
Q-balls of charge q, we have 

Q = Nq (3.24) 

and 

E = N{3 3 q 3/A = f3 3 Q 3/A N 1 ^. (3.25) 
Then energy considerations dictate 

/? 3 iV 1/4 Q 3/4 < PiQ 1/2 L, (3.26) 

i.e., 

(| ) V 2 < (^) 1/2 (3.27) 
which says the size of each small Q-ball must be greater than the width of the Q-wall. 



4. Stability 

The general rate of decay is governed by the stability analysis of the fluctuations 
about the Q-wall or Q-string configuration. The total Lagrangian is expanded to 
second order in the fluctuations about the soliton configuration. The first order 
terms vanish because the soliton is a solution of the full 3 dimensional equations 
of motion. The second order truncation yields a generalized normal mode problem. 
Real frequencies lead to oscillatory behaviour while imaginary frequencies yield ex- 
ponential growth or decay. In the latter case, the lifetime is taken to the inverse of 
the magnitude of the imaginary frequency, which is the time scale for the instability 
leading to beading. We will analyze the problem of the stability of a Q-wall for a 
specific potential which leads to an analytically accessible problem. 

4.1 Analytical Q- walls 

As we have seen in section |2.3| , the expression for a Q-ball in 1 dimension can be 
brought to quadrature. For many choices of the potential the resulting integral can 
be computed analytically 0. The stability analysis of the corresponding Q-wall 
is related to the elegant theory of exactly solvable quantum mechanical systems, 
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super symmetric quantum mechanics and shape invariant potentials ||. Consider 
the potential 

V{<p)=«?£<f?-2 ( £). (4.1) 
2 0o 

Then 

-V(<p) + = (4.2) 

implies 

^ = or cp = <p . (4.3) 

The non-trivial zero is fixed at 0o an d the Q-ball solution will start from this value 
for (p at its centre. The equation of motion is 

-tff^x) - cu 2 <p (x) + V'(<p (x)) = (4.4) 

which becomes 

-<p%(x) + uo\A^{x) - 6^M) = 0. (4.5) 
The quadrature [2.19| is integrable in this case, with solution 



(po(x) = 0osech (ux). (4.6) 

It is remarkable that the stability problem around such a Q-wall is related to the 
corresponding suspersymmetric quantum mechanical system, as we will now show. 

4.2 Second-order Lagrangian 

The stability analysis proceeds with the second-order expansion of the total La- 
grangian about the Q-wall configuration. We will consider an infinite planar Q-wall, 
which then gives that ^4.6| is an exact solution of the full 3 dimensional equations. 
This will be an arbitrarily close approximation for a large but finite closed Q-wall sur- 
face. The second-order Lagrangian for 4>(x,t) = e luJt (p(x,t) = e luJt (<p (x) +<p\(x,t) + 
iip 2 {x,t)) is 



/l 11 
d 3 x-(d t (p a -uje ab (p b )(dt(p a -uje ab ip b ) + -di<p a di<p a - -V"(ip )(pl - 



l%o) 2 
2 if ^ 2 
(4.7) 

where a, b = 1, 2. Arbitrary variations <pi and <y? 2 are permitted from the point of view 
of charge conservation. All spatial Fourier components in the orthogonal directions 
will correspond to displacement of charge along the wall in a periodic fashion. Charge 
conservation constraints will only apply to the zero frequency Fourier component, 
that is constant along the wall, and here the variations will have to be constrained. 
But this component corresponds to fluctuations that are effectively 1 dimensional, 
depending only on the coordinate normal to the Q-wall configuration. We know 
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that the Q-wall is by definition stable under such fluctuations, hence we will not be 
concerned with these fluctuations. The equations of motion for ip a are 



<fl - 2uJ(p 2 - C^Vl _ VVl + V"((f )(f 1 = 

(p 2 + 2u<p x - wV 2 - VV2 + V — <f2 = 0. (4.8) 

Replacing ip a (x,t) — > e~ % ^ nt+klV+k2Z ^ (p a {x) and writing the resulting equations in ma- 
trix form yields 



. l2M ^ t -M. ww =f-«g («) 

which is a transcendental equation for f2. 

The operator on the LHS plays the role of a Hamiltonian whose eigenvalues 
determine if a mode is oscillatory or exponential. The Hamiltonian is a hermitean 
operator for real Q, hence its spectrum is real. It is easy to see that all non-negative 
eigenvalues of the Hamiltonian at Q = will give rise to oscillatory modes of the 
fluctuation spectrum. Only negative modes at O = 0, for small enough values of 
I A; I can yield exponential modes, but even these become oscillatory modes for large 
enough \k\. The non-standard aspect of this stability problem is that the Hamiltonian 
also directly depends upon the fluctuation frequency Q. We must find the spectral 
curves E(Q) as functions of Q and then look for solutions of the equation E(Q) = 
Q 2 — \k\ 2 . The RHS is an upward opening parabola for real Q with intercept at 
— |A;| 2 . The spectral curve of any eigenvalue of the Hamiltonian that is non- negative 
at Q = will intersect with the parabola at two distinct points. This is because 
for large enough Q the spectral curves are simply asymptotic to the linear functions 
E(Q) = ±2u;Q. It is only the eigenvalues of the Hamiltonian that are negative at 
Q = which can yield exponentially decaying solutions. We will see that there is 
exactly one such mode for this Hamiltonian with the choice of the scalar potential 
taken in 14.11. 



4.3 Spectrum of the Hamiltonian and supersymmetric quantum mechan- 
ics 

The Hamiltonian for Q = is given by 
There is an obvious exact zero energy level, 

.<P2/ \Vo(x 



r = rVf..!.,] (4.11: 
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since (fo satisfies the equation |4l|. Remarkably the mode 



^.V.WrtM, (4 . 12) 



(f2 ) V 

is also an exact zero mode, essentially the corresponding eigenvalue equation is the 
derivative of the equation of motion for (fo(x) given by equation [4.5| . The nota- 
tion will become clear below. This is, of course, not just a coincidence. The two 
sub-Hamiltonians which comprise Ti are actually supersymmetric partners and their 
spectra are paired except for the lowest energy mode. Indeed since |4.12| is the deriva- 
tive of ipo(x) it has a node, hence it cannot be the lowest energy eigenvalue of TL + , 
and there has to be a negative energy mode. 

4.3.1 Stability equations and supersymmetric quantum mechanics ||] 

Consider the class of supersymmetric quantum mechanics models determined by the 
super potential 

W = nujtanh(ujx). (4.13) 
Then the corresponding supersymmetric partner Hamiltonians are given by 

H± = W 2 =F W = nV - n(n ± n)u 2 sech 2 ux (4.14) 

The two Hamiltonians 7i± are given in terms of 

H+ = A ] A H- = AA + . (4.15) 



where 



which gives 



A = ^- + W A^ = -^ + W (4.16) 

dx dx 



d 2 

H.+ = — + n 2 ^ 2 — n(n + l)d; 2 sech 2 a;x 
dx 2 
d 2 

H- = — + n2 ^ 2 - n(n - l)u 2 sech. 2 ujx. (4.17) 
dx 2 

The ground state is the wave function annihilated by A 

Atp-x = (4.18) 

which integrates easily to 

(p-i(x) = </?_i(0)sech T W. (4.19) 

Then, if 7i + ip + = E + ip + then 7i-{Aip + ) = E + (Aip + ) which shows that all energy 
levels are (supersymmetrically) paired except for the level annihilated by A. Equiv- 
alently if 7^_^- = E_ip_ the 7i + (A^il)_) = EL(j4'^>_). This pairing is equivalent to 
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the former except that there is no state annihilated by A^. To find all the levels of 
7i± we proceed by changing the value of n — > n — 1. This process continues until we 
reach n = 1 and the Hamiltonian 7Y_| n =i = — (d 2 /dx 2 ) + uo 2 for which we know the 
spectrum to be free, plane waves. The ground states of all intermediate Hamiltonians 
are given by 

ip ~ sech"W m=l, 2,---,n. (4.20) 

Then applying the appropriate strings of A]^ (using an obvious notation) generates 
all of the bound states of TC± as shown in detail in ||. The interesting point is that 
we can easily find the potential V(ipo) that is necessary to generate this system of 
super symmetric Hamiltonians for the stability problem of the Q-wall. Indeed, we 
impose that 

V '(fo) 2 / -n 1,2 ,r» 

= io nvn — ljsecn ux + B 

V"(if ) = u 2 n{n + l)sech 2 a;a; + C (4.21) 
where B, C are constants to be adjusted. This yields 

1 f V (<p ) 5) = w 2 nsech2LJX = _J_(yfi^ _ C y ( 4> 22) 



n — 1 tpo n + 1 

This equation is equivalent to 

(4.23) 

after multiplying by an integrating factor. Then this equation is easily integrated to 



VM = Rg^ + ( ^) : C(n-l) H + G (424) 



for some integration constants D, G. Now imposing consistency with equations |4.21 
yields 

Y1M. = Dip ^ + ( B(~ + D-c( n -i) = w2n(n _ 1)sech2aii + B 

Lf 2 

Taking = u 2 n(n — 1) yields 

Vo = sedi^ux and B i n + l ) ~ C{n-\) = R 



On the other hand from |4.21 



,r,u \ + -( ZLt r)- 1 /-B(n + 1) - C(n - 1), 2 , „ N l2 ^ 
^"(^o) = D(— yH + (— ^2 — " -) = w n (™ + l)sech 2 wx + C 

(4.27) 



15 



hence taking D = uj 2 n{n — 1) as before is consistent since 

n 4- 1 

£)(____) = ^ n ( n + i) (4.28) 
n — 1 

as required. The constants however must satisfy (along with the equation from 
before) 

B{n + 1) - C(n - 1) D 

g <" +1 )^ C '"- 1 ' = C. (4.29) 

These are two homogeneous linear equations for B, C which only have the trivial 
solution except if the system is dependent. It is indeed, and we find the solution is 
simply B = C. We see that the potential is a polynomial for only the cases n — 2, 3. 

4.4 Stability analysis for analytical Q-walls 

The system studied in section [O] uses 

V'{ip {x)) = ( - 6u 2 ^P- + 5u 2 )Mx) (4-30) 

hence replacing <po(x) = 0osech 2 o;x yields 

YM*» = (5. 2 - 6,/^) = ^ - 6^sech 2 ^ (4.31) 
(p {x) 0o 

while 

V"((p (x)) = {5u 2 - 12o> 2 ^^) = 5o> 2 - 12cu 2 sech 2 cux. (4.32) 

00 

Thus the system corresponds to the case n = 3 above with B = C = 5o> 2 . The 
corresponding stability equation is 

—£s + 4a; 2 - 12o; 2 sech 2 a;a; i2ujVt \ ( <p x 

-i2ujQ --£; + 4a; 2 - 6o; 2 sech 2 a;a; J \<f2 

(n 2 -\k\ 2 ) . (4.33) 

The spectrum of this Hamiltonian at Q = corresponds to a single negative energy 
bound state with E = —5a; 2 



15a; / sech uix 
16~ V 



^ = \l—[ n ) ( 4 - 34 ) 



two zero energy bound states with E = 



!l5u ( sech 2 o;xtanha;x\ , /3a; , 
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and two positive energy bound states with E = 3uj' 

.3 



+ prr— ( 12sechu;x — 15sech ux \ 3u ( \ , . 

^ 3uj2 ~ U \ J ~ V~2~ Vsech^xtanhcux / 

and a continuum starting at E > 4uj 2 . All these modes are found in a straightforward 



manner using the general ideas schematized in section [4.3 . 1|, for more details see 



The eigenvalue problem [4.33| has solutions that are best expressed graphically. 
The spectral curve E(Q) of each eigenvalue of the Hamiltonian on the LHS of |4.33| 
must intersect the parabola Q 2 — \k\ 2 . The eigenvalue problem for this Hamiltonian 
can be written algebraically as 

(H + #i<t 3 + 2uQa 2 )ip = E(Q)if) (4.37) 

where <7j are the usual Pauli matrices. The spectrum is doubly degenerate for each 
non-negative mode at Q = due to the underlying supersymmetry. Clearly 

(Ho + H i( t 3 - 2uVL<j 2 )a^ = E(n)a 3 ip (4.38) 

hence the spectrum is even in Q. However we can without loss of generality restrict 
Q > because <fi are in principle real: reversing the sign of Q — > —Q does not 
generate an independent set of solutions. 

We see clearly that the non-negative modes satisfy E(Q) > Q 2 — \k\ 2 for Q 0, 
but as Q — » oo these spectral curves asymptote to E(Q) — > ±2uQ, the spectrum of 
the operator 

2uQa 2 . (4.39) 

These are linear functions of Q which for sufficiently large Q always satisfy Q 2 — \ k\ 2 > 
2uQ. Hence by continuity each spectral curve must intersect at least once with the 
parabola, giving rise to a real solution for Q, and hence an oscillatory mode in the 
fluctuation spectrum. We speculate that this occurs exactly once for the spectral 
curve of each non-negative eigenvalue. 

Only the negative mode needs to be treated differently. Here, at Q = and for 
\k\ 2 small enough, the negative mode is already below the parabola and hence there 
is generically no real solution for Q. We cannot compute the solution analytically, 
but a perturbative approach is fruitful. We can compute the spectral curve to second 
order in perturbation theory (the first order change vanishes). We get 

£l(n) = -5o, -40 u; ^ + ••• 

n=l 

= -5u 2 - Q 2 B 2 + o(0 4 ) (4.40) 

where B 2 is a calculable positive real number that is independent of uj. It is not 
illuminating to compute B 2 , it is a number of order 1. Then the equation to be 
satisfied is 

-5uj 2 -n 2 B 2 = n 2 - \k\ 2 (4.41) 
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0.2 0.4 



Figure 1: Graphical depiction of transcendental equation for ft (Eq. [4. 33 , with an appro- 
priate numerical choice for oj); k 2 > .5. Solid line represents discrete eigenvalues of the LHS 
of this equation; continuous eigenvalues are expected to have similar behavior. Broken line 
represents the RHS. Points of intersection give the (real) values of O for the corresponding 
modes; no intersection would indicate no real solution, indicating an unstable mode of the 
Q-wall. No such instabilities exist in this case. 




0.2 0.4 



Figure 2: Graphical depiction of transcendental equation for Q, (Eq. [4. 33 , with an appro- 
priate numerical choice for to); k 2 < .5. No solution exists for lowest eigenvalue of the LHS, 
indicating an unstable mode. 



that is 



\k\< 



For | A; | 2 < 5uj 2 this gives Q 2 < 0, i.e., the solution for fl is imaginary and the 
corresponding mode is unstable. Indeed 



n 



5- 



%IQ\ 



l + B 2 



(4.43) 
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hence as \k\ 2 — > 5uj 2 from above, we get Q — > 0. The lifetime is defined as T = l/|f2| 
hence for small values of \k\ we get T ~ (1/^)- But as \k\ 2 — > 5c<j 2 we get T — > oo. 
Thus the very long wavelength spatial Fourier modes decay with lifetime of the order 

— * 

of 1/u while as the wavelength l/\k\ approaches l/5u, which is of the order of the 
thickness of the Q-wall, the lifetime becomes arbitrarily long, the modes become quasi 

— # 

stable. For larger \k\ (shorter wavelength) the modes are of course strictly stable. 
Analysis of related, exactly solvable simple models such as H{VL)ip = {k 2 + acr 3 + 
Q<j2)ip = Q 2 ip yields the same kind of behavior, generically. Hence we can conclude 
that the fluctuation spectrum corresponds to stable oscillatory modes except for one 
degree of freedom. The corresponding mode is unstable only if the wavelength of 
the spatial Fourier mode is larger than the width of the Q-wall. Furthermore, the 
lifetime of the unstable mode is larger than \/uj ~ Q^ 2 ^ = (Q / L 2 )^/ 2 \ which clearly 
can be taken to be arbitrarily large. Hence we conclude that Q-wall, or equally well 
a Q-string can be arbitrarily long-lived. 



5. Discussion and conclusion 

In this paper we have studied lower- dimensional Q-balls and extended objects based 
on them in models with a scalar field potential with a flat direction for large field 
values. The energy of a Q-ball of charge Q d in d space dimensions varies as Q^ d+1 ^ . 
Thus, as three-dimensional Q-balls become more and more energetically advanta- 
geous with increasing charge, the same is true (only more so) in two or in one 
dimension. 

The extended configurations, called Q-strings and Q-walls, also have interesting 
energetics. A string of fixed charge and of varying length L, for example, has energy 
E ~ L 1 / 3 , indicating a rapidly decreasing tension as the length is increased. The 
time scale of oscillations turns out to be of the same order as that for a topological 
string. 

These strings and walls have a different instability not shared by topological 
strings: namely, density fluctuations along the string or wall can form, causing a sort 
of beading: the extended object would evolve into a one- or two-dimensional array 
of smaller Q-balls. 

This instability was examined in detail for Q-walls, where supersymmetric quan- 
tum mechanics enables a rather detailed analysis. Only modes corresponding to 
long-wavelength fluctuations are unstable, and the time scale of growth of these 
fluctuations can be rather large (compared to the characteristic time scale of the 
Q-wall) . 

It is not clear whether or why such objects would form in the early Universe. 
(Unlike conventional cosmic strings, there is no topological impetus for their forma- 
tion.) But if they did form, they could be long-lived, and could help contribute to 
structure formation in the Universe. Such a possibility merits further investigation. 
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